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Abstract. Evolution with mass segregation and the evolution of the 
rotation of cores are both discussed for self-similar core collapse. Evolu- 

1/2 

tion with VLq oc pJ is predicted. On the Dynamical Main Sequence of 
globular clusters the energy emission from binaries balances the energy 
expended in expanding the halo. 

Newton's exactly solved N-body problem is then given, along with 
recent generalisations, all of which have no violent relaxation, but a new 
^ • type of statistical equilibrium is discussed. 

! Finally, we set the creation of streams in the Galaxy's halo in the 

historical context of their discovery. 

o 
o 
o 
o 

■ 1. Introduction 

When, over ten years ago, the main focus of my work moved away from stellar 
» ■ dynamics, three outstanding problems were left unsolved and my efforts to in- 

terest students in these have not borne fruit. In the hope that others may be 
stimulated these questions are posed here. 

Ql. Is it possible to have self-similar core collapse when there is a continuous 
^ ■ distribution of stellar masses? 

Q2. For a weakly rotating cluster is there a self-similar collapse of the core 
and how does its rotation evolve as the core radius, r c , decreases? 

Of course question 2 can be tackled for equal mass stars and then one can 
make it more realistic by combining both questions together in Ql + Q2. 

Q3. It is now widely accepted that just as in stars the energy generated 
by nuclear reactions delays thermal evolution, so in star clusters the energy 
generated by binary stars delays core-collapse. Eddington's theory of the Main 
Sequence led to many observational consequences. Can we make a theory of the 
Dynamical Main Sequence of globular clusters with real predictive power as to 
the form these clusters should take at equilibrium? 

Section 2 outlines possible solutions but full solutions involve multidimen- 
sional problems in four or five dimensions (e.g., Energy, Angular Momentum, 
Time and Mass). 

Some years ago Inagaki and I (1990) showed that such problems could be 
tackled analytically using trial functions and a Local Variational Principle. I still 
think that method should have a great future in these problems but those wishing 
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to use it should contact Takahashi (1992,3,5,6) who developed this approach 
numerically for 3 dimensional problems, but has since turned to another method. 

Since few workers in N-body dynamics have even heard of Newton's exact 
solution of an N-body problem for all N and all initial conditions, it is given in 
section 3 along with recent developments of this idea. These problems lead to 
some pretty dynamics which illustrate that there are exceptions with no violent 
relaxation and that thermodynamic equilibria can exist in a system undergoing 
a macroscopic radial pulsation or expansion. 

Finally in section 4 we discuss the history of the discovery of tidal streams 
and give the methods recently invented for finding new ones. Once some globular 
cluster extensions or proper motions are reliably determined these methods will 
become far more powerful. 

2. Core Collapse, Mass-Segregation, Rotation and the Dynamical 
Main Sequence 

Much of globular cluster evolution follows from careful dimensional analysis, 
see Inagaki & Lynden-Bell (1983). Returning to our Question 1 we look for a 
self-similar solution with the density in the core and inner halo of the form 

p(r, t) = po(t)p* (r*) where r* = r/r c (t) 

Here r c (t) is the core's radius. 

Since the relaxation is far more rapid in the core than in the outer parts we 
may set d/dt ~ for large r. 

= (dinp/dt) r = d£npo/dt — (d£np*/dr*) dtnr c jdt 

Therefore for large r 

d£npo/d£nr c = dlnp*j dtnr* = —a 

Since the expression on the left is a function of t alone, while that on the 
right is a function of r* alone, it follows that both are constant (at a value we 
have called —a). Thus for all time 

po = Ar~ a (1) 

and at large r*, = r~ a . 

The evolution of the core is due to the relaxation so 

Mnpo/dt = K 1 /T c (2) 

where K\ is a dimensionless constant and T c is the relaxation time in the core. 
A standard evaluation of T c gives for a velocity dispersion a, Binney & Tremaine 
(1987) 

T' 1 = 3G 2 m P en A /a 3 (3) 

and we would like to evaluate this for the core when there is a distribution of 
stellar masses. It would seem natural to take the distribution function to be an 
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equilibrium Maxwell-Boltzmann there proportional to exp —(5m {^v 2 — tpo^j 

but this cannot be true at energies close to the escape energy. The closeness of 
the escape energy is well emphasised by the following argument from the Virial 
Theorem. 

2T = Y,m lV j = -V = \Y, m ^i ( 4 ) 
I 1 I 

Now ipj = 7}V 2 where v e is the velocity of escape to infinity, hence for any cluster 
obeying the virial theorem 

(mv 2 ) = |H 2 ) . (5) 

So that the mass weighted rms velocity is HALF the mass weighted rms velocity 
of escape. At a thermodynamic equilibrium we have equipartition of \mv 2 so 
the above expression implies that stars with less than a quarter of the mean 
mass have an rms velocity equal to the escape velocity. This cannot occur in 
practice. Thus the equipartition of the Maxwell-Boltzmann distribution must 
be modified and as Michie (1963) first showed the lowered Maxwellian of the 
Michie-King models [evaluated in full by King (1966)] does this approximately 
viz 

-(3m I — v 2 — ip 



f = B{m) < exp 



2 



exp {(3m%l) e ) \ (6) 



Although in well developed core collapse such models have insufficient temper- 
ature contrast between the centre and the halo, Lynden-Bell & Wood (1968), 
Lynden-Bell & Eggleton (1980), nevertheless we shall adopt such a form for the 
low energy stars in the core (for which relaxation is the most rapid). If F(m)dm 
gives the number density of stars at mass m in the core, so that 

po = J mF(m)dm (7) 

then we may re-express the relaxation rate in terms of the central 'temperature' 
(3~ l = o 2 jm and obtain from (3) 

T~ l = 3G 2 l3 3/2 p £n A (m 7 / 2 )/{m) (8) 

where 

V 2 ) = J m 7 / 2 F(m)dm/ J F{m)dm (9) 

Evidently the relaxation rate, T" 1 , which determines the overall evolution 
of the cluster, depends on the mean seven halves power of the mass, evaluated 
over the core of the cluster. As evolution proceeds, the lighter stars in the 
core are preferentially expelled from it but simultaneously the heavy stars are 
gradually eliminated via stellar evolution. Before using (8) in (2) we need to 
re-express (3 in terms of po, r c and (m). If M c = ^irpor 3 is the core mass then 

3a 2 = 3((3(m)y l = GM c /r c = —Gp r 2 (10) 



(m 
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where we have chosen to define r c so that the constant of proportion is one. 

T- 1 = K 2 G^Po 1/2 r- 3 m c = K 2 G X I 2 A^m^l Af>-<*IV<*) (11) 
where we have used (1) to express r c in terms of po, K 2 is the dimensionless 
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-£nA) and m c (t) = (m 7 / 2 ) / (m) 5 / 2 . To use expression (11) to 



constant , 

V Sir- 1 ' * 

solve (2) we must first evaluate the time dependence of m c . This is caused by 
two separate effects, the progressive expulsion of lower mass stars from the core 
increases m c but it is mainly dependent on the high masses and they diminish 
steadily as stellar evolution and stellar death take their toll. With such a high 
power of the mass involved it is likely that stellar evolution plays the dominant 
role with m c behaving similarly to the age cut off. This suggests the approximate 
form 

m c = m Q (t/to) 5 ' 1 where S ~ 2/3, t Q ~ 10 10 yr, t > 10 6 yr 

and t is time measured from the time the cluster was created. Inserting this into 
(11) and (11) into (2) the resulting equation for po is readily solved to give 

/9o = C|4-t 5 j- 2 «/(6-a) (12) 

and therefore 

r c = (A/C) 1/a \t s c -t s \ 2 ' {% -^ (13) 

where t c is a constant of integration which gives the time of core collapse and 
the constant C is 



C = A 



^K^A-^m^ 



-2a/(6-a) 



As in Inagaki & Lynden-Bell (1983), we expect (12) and (13) to hold also after 
core collapse but then po becomes a characteristic density rather than the central 
one which remains infinite if we ignore binaries and the giant gravothermal 
oscillations. 

Of course in the absence of stellar evolution m c will increase and the 
parameterisation 

m c = m Q (t/T c ) s 

might then be more appropriate but though such a model is soluble I shall not 
detail it here. 

Expressions (12) and (13) can only be considered the solution to our problem 
once a is known. As Lynden-Bell and Eggleton showed a emerges as an eigen 
value in the full theory and for the equal mass case a = 2.22. It is readily seen 
that a must lie between the equal mass isothermal sphere with a = 2 and the 
limiting case with infinite core binding energy a = 5/2. With the heavy masses 
more concentrated to the core one may expect 2.22 < a < 2.5 and perhaps 
toward the upper end of that range, on the other hand once most of the lighter 
stars have been ejected from the core one expects the evolution to revert toward 
the single mass case. Our best guess is therefore 



a = 2.3 ±0.08 



(14) 
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In discussing Question 2 we note that the flattening of a cluster is of second 
order in Q the angular velocity and may be neglected to first order. In any 
potential of the form 

i/> = ik(r)+r- 2 M0) ( 15 ) 

the expression / = ^ J 2 — ^2 is an exact integral of the motion where J = (r x v) . 
In practice any non-spherical part of ip is second order in f2 and the motion of a 
star is well approximated as lying in a precessing plane. The rate of precession 
can be worked out from the couple that the non-spherical potential exerts on the 
unperturbed rosette orbit that the star would have in the absence of asphericity. 
Thus in practice | J| is almost an integral of the motion while J z and the energy, 
e, are exact ones in the absence of evolution. In practice it is convenient to use 
the radial action J r = 4- § \J1\e + ipo(r)] ~ J' 2 /r' 2 dr in place of e. In potentials 
of the form (15) it is an exact integral but for more general forms of ip it is 
only approximate. A particular advantage of using |J|, J r , J z is that they are 
adiabatically invariant for slow changes in the potential ip. Thus if we express 
the distribution function of a rotating cluster in the form / = /(|J|, J r , J z , t) 
where / evolves little in one orbital time, then df /dt is (df/dt) encounters and 
there is no supplementary term due to the change of global potential induced by 
these encounters since for such changes the J (although not e) are adiabatically 
invariant. The J have dimensions of or c and / has the dimensions of po<J~ 3 , so 
in self-similar evolution the evolving / will have the form 

f = p a- 3 F*(\r\, J r *, J*) (16) 

where is a dimensionless function of its dimensionless arguments 

| J* | = |J|/(crr c ), J* = J r /(ar c ), J* = J z /(ar c ). As before p = Ar~ a so 

or c oc r 2 a ^ 2 ■ These same principles may be applied to discuss the way cluster 
core rotation evolves as the core contracts. Near the centre, relaxation is quite 
fast so we expect a rotating Maxwellian at energies little affected by escape 

f = Bexp-a- 2 (e-Q J z ) . (17) 

Our question is how does f2o evolve during core collapse. In self-similar evolution 
the quantities that are dimensionless do not evolve, so £lor c /a = const. Hence 
as the core shrinks we have our prime result 

ft oc r- a/2 a ~ 2.22 . (18) 

So as r c shrinks, the core should rotate faster and faster. This result is very 
sensibly between the ^o r c = const that would follow if the heat flowed out of 
the core much more readily than angular momentum, and the 0, = const that 
would follow if the opposite were true. Away from the centre it is plausible 
that the system remembers the rotation it had when that part of the inner halo 
was part of the core; in reality it gains a little more angular momentum as it 
leaves the remaining core, because the remaining core loses it. This suggests 
the behaviour ~ r~ Q / 2 but since VL 2 and Gp have the same dimensions a 
better prediction is that behaves approximately as p 1//2 ; hence we conclude 
that cluster cores and inner haloes rotate as 



SI = Sl (p/po) 1/2 = QoP* 



(19) 
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with Qq evolving according to (18). For Q\ + Q2 the same results will hold 
with a given by (14). Such predictions should be compared with the numer- 
ical work of Spurzem (2000) and the rapidly growing data from observations. 
A further consequence of (19) is that SI 2 /itGp which determines the flattening 
will be constant in the core and inner halo thus the isophotes there should have 
b/a constant whereas if the core evolved with CIq oc r~ 2 the central isophotes 
would be more flattened. By contrast a uniformly rotating cluster has the great- 
est flattening in the outer isophotes. Turning now to Question 3 the theory of 
the Dynamical Main Sequence of Globular Clusters as yet awaits someone brave 
enough to make strong hypotheses such as 'The only binaries that matter are 
in the core and the energy flux through the inner halo is constant and drives 
either the escape from the cluster or the expansion of its halo.' Of course if 
core collapse ceases due to binaries then the steady core will relax to uniform 
rotation which will gradually extend from the core outwards. 



3. Newton's N-body Problem and its Generalisations 

Let the force on body I due to body J be of the form 

F/j = fcm/mj(rj - 17) (20) 
We sum over all J to get the force on / since the I = J term is zero. 

F/ = kmj(^2 m J r J — Mrj) = kmjM(f — 17) 

where M is the total mass. Thus the total force is directed to the centre of mass 
r and is proportional to the distance from it. Newton's (1687) general solution 
is that each particle moves in a central ellipse about the centre of mass which 
itself moves uniformly in a straight line. Newton's system has a total potential 
energy 

V + -KM 2 r 2 (21) 



where 

The Virial theorem reads 



r 2 = M-^mKrj - f) 2 (22) 

-I = 2T + nV (23) 

where V oc r~ n so n = — 2 for the above system. To generalise. Newton's result 
consider systems with total potential energies of the form 

V = V{r) 

with r given by (22). This problem has a singular beauty if we consider the 3N 

dimensional space in which the X/ = ^J^jfi^i — f ) are the coordinates. Letting 

the 3N vector r = (xi,X2 • • • , xjv) we note that r 2 is given by (22). The initial 
conditions in this space are the initial values of r and r. The accelerations in 
this space are central since V = V(r) so the acceleration also lies in the plane 
defined by the initial r and f . Thus the motion continues in that plane. In fact 
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the whole problem now reduces to the planar orbit problem under the action 
of the central potential V{r). Back in 3 space each particle feels the central 
force V'(r)miM~ 1 (rj — r)/r so each particle orbits in a plane. Its motion may 
be obtained by projecting the motion of the representative point in 3N down 
into the plane of the motion of particle I. As r changes periodically around the 
planar orbit in 3N space it follows that Mr 2 vibrates forever so there is no violent 
relaxation to a Virial equilibrium. One may generalise this result still further 
by taking V = Vb(r) + r^ 2 V 2 {Y/r) then the Virial theorem reads 



-Md 2 r 2 /dt 2 



2T - rVo(r) + 2V 2 = 2E - 2V - tVq 



since this last expression does not involve V2 it follows that r vibrates in the same 
way as it did when V2 was 0. Hence once again there is no violent relaxation. 
Pretty N-body problems of this type are given by the forces 



Fjj = Gm I m J {Y J - 17) 
where r is given by (22) or alternatively 



1 



A: 



(rj - rjY 



\-4 



Fu = G'mimj{r,j - 17) 1 - k{rj - 17) 

both of these force laws have long range attractions and short range repulsions 
so they may produce liquid-like and solid-like phases . In each case the vibration 
in r separates from the other motions so if lattices are possible, lattices with a 
breathing pulsation will be too. It is possible to do a statistical mechanics of 
the 'angular' motions only while r remains breathing but for details the reader 
should refer to Lynden-Bell &: Lynden-Bell (1999a). The 1999b paper solves 
Schrodinger's equation for both spin Bosons and spin ^ Fermions giving the 
only known non-trivial three dimensional N-body solutions for interacting Bose 
or Fermi gases. The degenerate Fermi case gives a white-dwarf-like solution. 



4. Streams 1950-2000 

Streams about the Galaxy have much in common with the meteor streams left 
behind by comets in the solar system which were discovered earlier. However the 
galactic orbits are rosettes rather than ellipses and their planes precess about 
the galactic pole due to the aspherical potential. Berth Lindblad (1958) (1961) 
was one of the first to consider the theory of streams spreading from a common 
origin but already Eggen (1959) (1989) was on the march finding moving groups 
in the velocity space among both disc and halo stars. As data were refined 
Eggen changed his mind about the central velocities of some groups and this 
led to consequent changes in their membership which only increased scepticism 
among his critics. However Eggen was convinced he was onto something and 
claimed that it was much easier for others to criticise the inclusion of a star when 
new data showed its membership to be unlikely, than to discover the group to 
start with. 

The subject of streams really come to life with the discovery of the Magel- 
lanic Stream stretching more than 120° around the galaxy. But this discovery 
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was the culmination of a sequence. In 1965 Nan Dieter discovered a high velocity 
cloud at the South Galactic Pole. Soon afterward Hulsbosch & Raimond (1966) 
published their survey of such objects and showed that cloud to be strongly 
elongated. It was the technological innovation of Wrixon's low noise diodes that 
allowed Wannier & Wrixon (1972) to see fainter 21 cm emission. Their great arc 
of 21 cm emission through Dieter's cloud showed a systematic variation of radial 
velocity by several hundred km/sec across the sky. At Herstmonceux Kalnajs 
(1972) excitedly drew my attention to their paper claiming that it must be of 
Galactic scale. A day or two later he discovered that both the line of the arc and 
the radial velocity extrapolated to the Magellanic Clouds. Mathewson too was 
very excited by this discovery and no sooner had he got back to Australia than 
he, Cleary & Murray (1974) started delineating how the stream joined the Mag- 
ellanic Clouds. A few years later during a non-photometric night at Sutherland 
I was imagining how the Magellanic Stream lay across the sky when it struck 
me that other satellites of the Milky Way might be members of streams and 
might likewise have high velocity clouds associated with them. I rapidly found 
that Draco & Ursa Minor are almost antipodal to the LMC and the SMC in the 
galactocentric sky. This meant that they too might be long lost members of the 
Magellanic Stream, Lynden-Bell (1976). Kunkel & Demers (1977) noticed the 
dwarf spheroidals were distributed non-randomly over the sky but their com- 
promise plane was some 30° from the Magellanic Stream. In their attempt to 
explain the antipody of Draco & Ursa Minor, Hunter Sz Tremaine (1977) noticed 
that their elongations lay along the stream while I searched for further streams 
using elongations as a supplementary guide, Lynden-Bell (1982); however it was 
not until 1995 that this developed into a systematic method. Since tidal tearing 
occurs essentially in the plane of orbit, one needs to find objects scattered along 
a great circle in the galactocentric sky. Such great circles are hard to see in 
any of the projections of the sky so we need a good method of finding them. 
An object may be associated with the set of all great circles that pass through 
it. The poles of these great circles lie on another great circle which we call a 
polar path. (The pole of that great circle lies of course at the object we started 
with). The great circle through any two objects has its pole at the intersection 
of their polar paths. If 3 objects lie on a great circle all 3 polar paths will go 
through the same point. Thus we may find objects with common great circles by 
plotting their polar paths (in any equal area projection) and looking for multiple 
intersections. Figure la which plots 45 polar paths for dwarf spheroidals and 
outer globular clusters shows a result that more closely resembles a ball of wool 
than a useful scientific plot. More information is needed to get any definitive 
groupings. A known proper motion would reduce the polar path of an object to 
a very short segment determined by the directional error of the proper motion. 
If we are prepared to assume that the elongation of an object is along its orbit, 
as appears to be true for Draco & Ursa Minor and might be true for the outer 
extensions recently found around globulars (e.g., Grillmair et. al., Meylan this 
volume) then likewise the pole must lie in only a small sector of the polar path 
(see Figure lb). But great circles are not enough. Radial velocities provide a 
discriminant. If all the objects in a stream follow the orbit of the progenitor 
then they will share the same specific energy E and specific angular momentum 
h. 
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45 < r g < 151 kpc 



Equal area projection of the sky b > -12deg. 
Polar paths of all 45 objects do not by themselves 
show obvious streams. 

X = Vl-sin bg cos tg, Y = Vl-sin b g sin lg where 
(£g,bg) are galactocentric galactic coordinates. 
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As Figure la but with only distant objects 
plotted and those portions of the polar paths 
perpendicular to the major axes of the 
objects emphasised. 



Figure 1. a & b 




Figure 2. The radial energy plotted against r g 2 . The gradient fives 
—h 2 /2 where h is the specific angular momentum of the stream. From 
this the predicted proper motions are deduced (see L-B) 2 1995. 
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Taking sets of objects that share a common polar intersection and therefore lie 
on a galactocentric great circle one may plot E r = E — ^h 2 r~ 2 = ^v 2 — ip{ r ) 
against r~ 2 the inverse square of their galactocentric distances. This should be 
a straight line so E and h 2 can be deduced as the intercept and gradient of such 
a line. In practice the observed radial velocity is not the galactocentric one so v r 
can be deduced by estimating an initial h and then iterating to find a consistent 
gradient of — \h 2 . In practice this method converges rapidly but needless to say 
some of the objects do not lie on straight lines and are interlopers rather than 
stream members. Since the method gives h? from the observed radial velocities, 
and since the poles of the great circle involved are known, we have used these 
methods not only to produce possible streams among the dwarf spheroidals and 
outer globulars but also to predict their proper motions. The method works well 
for the Magellanic Stream and for the Sagittarius Stream which was discovered 
by Ibata, Gilmore & Irwin (1995). There is now great hope that we can identify 
orbital streams associated with the later merging events emphasised by Searle 
and Zinn. Further work may take us back to the orbits of objects that took 
part in the initial formation itself. The ELS (1962) picture had the bulk of 
the metal poor globulars formed in the fragmentation of the initial collapse, 
but following Ambartsumian we know that most of the stars would have been 
formed in looser aggregates and stellar associations which subsequently broke up 
and spread around the Galaxy. The idea that ELS pictured a perfectly smooth 
collapse is not in their paper but has been added by others who have not read it. 
Eggen (1959) had already worked on streams in the halo, the fragmentation of 
collapsing systems was treated by Hunter (1962), and I was working on the Large 
Scale Instabilities of shape that would make the collapse non-axis symmetric, 
Lynden-Bell (1964). The non- linear behaviour of these instabilities was later 
treated in the simple non-rotating case by Lin, Mestel & Shu (1965). Those who 
want to read for themselves what was thought at that time would do well to read 
my paper to IAU Symposium No. 31 on the Formation of the Galaxy, Lynden- 
Bell (1967). Streams from the initial collapse are less likely to be preserved than 
the younger streams from late additions, but finding correlated stellar motions 
from the initial creation of the Galaxy was what ELS was all about and their 
method of using abundances and adiabatic invariants as fossilised evidence on 
how the Galaxy formed is probably its most lasting legacy to astronomy. The 
merging of galaxies was suggested as important in the very percipient paper on 
bridges and streams by the Toomre's (1972) and White's (1976) demonstration 
of fragmentation of a uniform overdensity showed that even within the free 
fall time of the whole much fragmentation and subsequent merging of those 
fragments occurs. 

Returning to streams there is too much modern work to detail all of it, but 
work by Johnston (1998) holds out the prospect of more accurate determination 
of the Galaxy's potential using streams of stars just escaping from a satellite 
galaxy. By contrast, Helmi & White (1999) take such stars to be emitted with a 
not inconsiderable velocity dispersion so such streams would be harder to detect 
and would define the potential less accurately. It is not yet clear to me which 
view is more realistic. Lastly, Putman et. al. have at last detected a forward 
stream from the Magellanic Clouds, albeit somewhat skew of the main great 
circle. Currently the most plausible explanation is that it was gas around the 
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SMC that was torn but much of the forward stream ran into the LMC and the 
rest was diverted by it, however other interpretations are possible. 

Acknowledgments. The theory of rotation was developed after the confer- 
ence while writing up this contribution under the hospitality of IUCAA, Pune, 
India. 
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